Influence of Surface Waves on Plasma High Harmonic Generation 
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The influence of surface plasma waves (SPW) on high harmonic generation (HHG) from the 
interaction of intense lasers with overdense plasma is analyzed. It is shown, that the surface waves 
lead to the emission of harmonics away from the optical axis. These off- axis harmonics violate the 
parity selection rules found from ID models. Further, our investigations in the highly relativist ic 
regime point towards the existence of a new SPW generation process. 



High harmonics generated from solid density plasma 
surfaces may be about to become the tool that enables 
the next breakthrough in ultrafast science (TJ [2]. Ex- 
tremely intense XUV attosecond pulses are emitted from 
plasma surfaces, exploiting the non-linearity of relativis- 
tic physics. Because plasmas are used, the process is not 
subject to the limitations of the methods for generating 
attosecond pulses in gases [Ij. 

The generation scheme works as follows: an intense 
laser is incident on a solid surface, immediately ionizing 
the material. The resulting plasma reflects the light due 
to its overcritical density. If the field intensity is suffi- 
cient, the reflection process becomes strongly non-linear, 
yielding high harmonics together with the reflected light. 
Different mechanisms have been identified to be respon- 
sible for the high harmonic generation (HHG), some 
of them already apparent in the sub-relativistic regime, 
while others take over in the highly relativistic regime 
p-5j. All of these theoretical descriptions, including nu- 
merical simulations, were formulated in simplified one- 
dimensional (ID) models, assuming the incoming laser 
light as well as the reflected radiation to be plane waves 
[3-6J. Diffraction and focusing of the extremely broad- 
band harmonics radiation in real 3D space can be treated 
separately, provided the HHG happens "locally indepen- 
dently" at each point in the focal plane [7]. 

One of the first theoretical results about the surface 
harmonics have been selection rules concerning the par- 
ity and the polarization of the generated radiation [6j . In 
case of normally incident light, the rule can be derived 
as follows: The radiation source term is the transverse 
electron current in the skin layer, jy = pPy/j. Here, 
p = — e(ne — ne,o) is the charge density (ions are as- 
sumed immobile), Py is the transverse momentum com- 
ponent and 7 the relativistic 7- factor of the electrons. As 
the reflected radiation contains the same Fourier compo- 
nents as the source current jy, we now examine, which 
Fourier components the individual factors may contain. 
In the assumed ID geometry, due to the conservation of 
canonical momentum, we have Py oc Ay. Thus, to first 
order, it oscillates at the laser frequency and we may 
write Py = Pskin^"*'^^. Oscillations in the relativistic 7- 
factor and the charge density p on the other hand are 
determined by the relativistic ponderomotive force and 
therefore must have a period which corresponds to half 
of the laser period. They contain only even harmonics. 



which we denote hy p = p2uj and 7 = 72^;. Consequently, 
the product jy = ppy/j contains only odd harmonics, 
and so does the reflected radiation. In a similar way, 
selection rules can be derived for oblique incidence, a 
complete overview of this was given in Tab. 1 of Ref. [6j. 

The generation of surface plasma waves (SPWs) is a 
vital feature of laser interaction with overdense plasma. 
A SPW can propagate long-distances along the surface, 
while its amplitude falls off exponentially in a direction 
perpendicular to the surface [SUTT], The dispersion rela- 
tions of the SPW reads as [tH [12]: 



where ujp is the electron plasma frequency, cjgpw and 
/cgpw are the SPW frequency and wavenumber, respec- 
tively. It has been known for some time that a laser 
impinged on the vacuum overdense plasma surface can 
generate two counter-propagating SPWs, and this pro- 
cess of SPW excitation was termed the two-surface waves 
decay (TSWD) process [T3HT6] . Physically, the pondero- 
motive force of the laser pulse causes density oscillations 
at the surface. These density oscillations, at twice of the 
laser frequency, beat with the electron oscillatory veloc- 
ity, causing the laser pulse to produce a nonlinear cur- 
rent density which oscillates at one of the SPWs frequen- 
cies. For normal incidence of the laser, both SPWs os- 
cillate with the laser frequency itself, forming a standing 
wave. The SPW becomes parametrically unstable and 
energy flows into the SPW modes [16j. The expression 
for the growth rate of the SPWs excitation (cf. Eq. (26) 
in Ref. [16j) reveals that this process becomes resonant if 
one chooses initial plasma density to be Up = 2cJo, where 
ujp is the laser frequency. The resonance occurs due to 
the presence of the dielectric function €2^00 — ^ ~ ^p/^^o 
in the denominator of the growth rate expression. This 
resonant excitation of the SPWs in the laser-solid inter- 
action can significantly change the surface morphology. 
Thus, if present, the SPWs are expected to affect HHG 
greatly as it occurs from the plasma surface. However, 
to the best of our knowledge, the influence of SPWs on 
HHG has never been investigated so far. 

In this paper, we study the influence of SPWs on sur- 
face HHG on the basis of two sample PIC simulation 
runs. We thereby focus on normal incidence, as the in- 
fluence of the SPWs can be pointed out most clearly 
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Figure 1: Electron density profile (a) before and (b) after the 
SPW emerged. Simulation parameters are: ao = 1.7, ne,o = 
5nc. Frame (a) was recorded at simulation time ct = 4A, (b) 
at ct — 30A. 



here. The first siraulation is in the moderately relativist ic 
regime, where we have reason to believe that the analytic 
description of the TSWD is approximately valid. The 
second is in the highly relativistic regime, where HHG is 
supposed to be most efficient and we expect considerable 
deviation from linear SPW theory. 

Let us begin with the first example. For this simula- 
tion, we chose parameters similar to the ones in Ref. [14J. 
The laser is a plane wave, normally incident on the tar- 
get. Its temporal envelope is given by a linear ramp 
over three laser periods, up to a constant amplitude of 
ao = eA/mgC^ = 1.7. The simulation domain is 2D, 
i.e. all derivatives in z-direction are neglected. Further, 
the laser E-field was chosen to point in ^/-direction. The 
plasma density leaps from to ne,o = within one grid 
step at X = 0.5 A. Ions are fixed. 

Figure [T] shows the surface density profile at two dif- 
ferent times. Frame (a) was recorded at time ct = 4A, 
just after the laser had hit the surface. Here, the elec- 
tron surface, initially located at x = 0.5 A, has been 
pushed in uniformly by about a tenth of the wavelength. 
Strong variations along the surface are not observed ex- 
cept for some small noise, so at this point the simulation 
results can still be well described by a ID model. This 
changes later in the simulation, as we can see in frame 
(b), recorded 26 laser cycles later. Here, a SPW has 
emerged, spontaneously breaking the translational sym- 
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Figure 2: /c^ -resolved harmonics spectra, recorded at the left 
boundary of the simulation box. Simulation parameters are: 
ao = 1.7, ne,o = Sric- Frame (a) corresponds to the time span 
ct — ... 10 A, when the SPW has not yet emerged, and (b) 
to ct = 25 ... 35 A, when the SPW is fully grown. The colour 
scale corresponds to the logarithm of the spectral intensity, 
normalization is arbitrary, but consistent. 



metry. The wavenumber, measured from the data shown 
in the figure, is given by k^^^ ~ 1.1 /^o, very close to 
what was expected from the linear SPW dispersion rela- 
tion Eq. ([l]), /cspw = 1.15 /cq. 

Corresponding spectra of the reflected light are pre- 
sented in Fig. [2] Again, frame (a) corresponds to an 
early state, before the SPW builds up. It is seen, that all 
radiation is reflected at /c^y = 0, as expected from a plane 
wave reflected at a smooth surface. We also notice, that 
harmonics are generated at odd multiples of the laser fun- 
damental, in agreement with the selection rules explained 
above. Also, we see that all spectral lines are distinctly 
broadened compared to later times. This can in parts be 
understood due to the motion of the reflecting surface, 
especially during the time span, when the laser intensity 
is still increasing [17J. Frame (b) shows the reflected ra- 
diation at the state, where the SPW is fully developed. 
We see that it possesses a clear signature in the harmon- 
ics spectrum. Now, light is not only emitted in specular 
direction, but also at ky = Ik^pw^ where / is an integer, 
numbering the angular sideband. These wavenumbers 
correspond to angles 



= arcsin 



Ick 



spw 



(2) 



away from the optical axis, wherein m is the number of 
the harmonic. Thus, e.g. the 3rd harmonic is measured 
in specular direction, but also at angles of ^\ = 22° and 
1^3 = 47° away from the optical axis. 

Another remarkable observation is that although the 
specular light still consists of purely odd-numbered har- 
monics, the angular sidebands contain both even and 
odd harmonic numbers. This can be understood with 
the following perturbative approach. As in the ID 
case, the source current is determined by three factors: 
the charge density p, the transverse momentum py and 
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the relativistic factor 7. Assuming that the SPW in- 
duces a smah perturbation, we neglect its influence on 
7, writing 7 = 72^; as before. To the factors Py and 
p, we add a perturbation term representing the SPW: 
P = P2u; + ^^Pspw cos(A:spw^) 6"^^^ and Py = Pskine"'^^ + 
^Pspw cos(A:spw^) e~*^^ Forming the product, we see that 



jy = jy^cos{ksp^y) 



SPspw P2u 
72a; 



^Pspw Pskin 



-2iojt 



(3) 



wherein jy corresponds to the unperturbed ID case, thus 
representing odd harmonics in specular direction. The 
second summand represents the first angular sideband at 
ky = ±/Cspw Obviously, both even and odd harmonics 
are contained here. Off- axis even harmonics arise due to 
the density modulation associated with the SPW, odd 
ones due to the modulation of the momentum. To assess 
the higher angular sidebands, terms of the order of 5'^ 
would have to be taken into account, including the non- 
linearity of the SPW itself. 

Figure [2] contains another important detail. We notice 
that when the SPW has appeared, the intensity of the 
higher harmonic orders (m = 5, 7) is considerably de- 
pleted compared to the ID case. Thus, it looks as if the 
TSWD process competes with the HHG process, consum- 
ing parts of its energy. This may explain why the HHG 
with laser pulses in the duration range of 30 ... 100 fs has 
sometimes fallen behind expectations from ID models 
and simulations. For longer pulses, heating of the plasma 
suppresses the SPWs again. 

For the second example, we increase both the laser 
amplitude and, so as to avoid relativistic transparency, 
the plasma density. The laser amplitude is raised to 
ao = 7.2, the plasma density to ne,o = 36 nc- If we 
generously extend the resonance condition of the TSWD 
to allow for excitation at harmonics instead of only the 
laser fundamental, we write the dielectric term as €200^ — 
\ — uJ^I 4c(j^ , where ujra = ttiujo . For our chosen density 
this term vanishes at m = 3, causing a TSWD resonance 
at the corresponding frequency. Hence for these parame- 
ters, the dominant excitation of the SPWs could possibly 
happen at the third harmonic instead of the fundamental 
frequency. 

Simulation results are presented in Fig. [sj Frame (a) 
shows the absolute square of the electron density, Fourier 
transformed in ^/-direction. This way, SPW modes be- 
come immediately recognizable as clearly defined spots 
in the image. In this example we see that the fundamen- 
tal does not efficiently excite an SPW, as there is no spot 
around ky = k^. This can be expected as the growth rate 
for the cjspw = ^0 mode is expected to be very small. In- 
stead, a sharp peak is observed at ky^^ ~ 3.3 /cq- This 
value is reasonably close to what is expected from the 
non-relativistic dispersion relation for cu. 



spw 



3a;o, which 

is ky = 3.7 /cq. We therefore conclude that here the 3rd 
harmonic excites the SPW instead of the laser fundamen- 
tal. To the best of our knowledge, this excitation process 
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Figure 3: Results of the simulation in the highly relativistic 
parameter regime: ao = 7.2, ne,o = 36 ric- Frame (a) shows 
the absolute square of the Fourier transformed electron den- 
sity, recorded at ct = 18 A. The colour scale is linear, units 
are arbitrary. Frame (b) shows the /cy -resolved electromag- 
netic spectrum, recorded at the left boundary of the simula- 
tion box, during the time window ct = 10 . . . 25A. Here, the 
colour scale is logarithmic. 



is completely novel and has not yet been described in the 
literature. Further simulations carried out by the authors 
indicate, that this type of excitation process happens fre- 
quently in the relativistic regime and is not limited to ex- 
act hitting of the assumed resonance. The demonstrated 
sample however shows a particularly clear case. 

As we see from frame (b) of Fig. [sj this process also 
leaves its footprint in the /c^-resolved spectrum of the 
reflected field. Here, no sidebands at ky = ±/co or ky = 
±2ko are seen, but one at ky = ±3ko is clearly observed. 
As expected, this sideband again contains both even and 
odd harmonic orders. So e.g. the 4th harmonic order is 
measured at an angle of 1^4 ^ 56° from the optical axis, 
but not on the optical axis itself. 

Figure [4] shows the growth and damping of the SPW 
modes in both simulations, represented by an integral 
over the spatial density spectrum. Frame (a) corre- 
sponds to the moderately relativistic simulation. It can 
be readily seen that after a short period in the begin- 
ning {ct < 10 A), where both modes grow about equally, 
the SPW at the fundamental frequency (blue line) grows 
strongly, whereas other modes are damped. This is ex- 
pected as we had chosen the plasma density to satisfy 
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Figure 4: Growth and damping of SPW modes, measured as 
the spatial spectral intensity of density fluctuations, defined 
by |n|^ = 1/ n(x,^,t) e:>q){ikyy) dy\^ . The blue lines corre- 
sponds to modes in the range 0.9 < ky/ko < 1.3, for this 
/c-range the SPW frequency is roughly cjq; the red lines cor- 
respond to 3.0 < ky/ko < 4.0, so that cjspw ~ 3c<;o. (a) shows 
results of the moderately relativistic simulation (ao = 1.7, 
ne,o = Sric), (b) of the highly relativistic one (ao = 7.2, 
ne,o = 36 ric). The black dashed lines show the exponential 
fits, from which the SPW growth rate has been calculated. 



the plasma resonance condition for the generation of two 
counter-propagating SPWs close to the fundamental laser 
frequency. After the SPW has reached a certain am- 
plitude, the plasma is heated rapidly and the wave is 
damped. Results of the highly relativistic simulation 
are displayed in frame (b). Here, the situation is the 
other way round: Whereas the SPW at 3rd harmonic can 
grow quickly, the SPW at fundamental frequency does 
not reach a high amplitude at all. This can be under- 
stood, because with ne,o = 36nc, the density is detuned 
far from the resonance for the conventional TSWD at the 
laser fundamental (ng = 4nc). Consequently, the growth 



rate is small and thermal damping sets in, before it can 
reach a significant amplitude. On the other hand, the 
3rd harmonic is exactly at resonance, if one allows for 
the generous interpretation of the TSWD theory men- 
tioned above. Thus the SPW at this frequency can grow 
quickly. Note however, that relativistic effects may have 
a strong bearing on the TSWD process and the SPW 
dispersion relation. These qualitative discussions are ex- 
pected to be correct, however quantitative comparisons 
are difficult to be drawn using the non-relativistic theory. 
From the PIC simulation, we find the growth rate to be 
F ^ 0.05 ujQ within the time range indicated in Fig.|4j Be- 
fore this stage of linear growth, there is a stage of quasi 
immediate growth, happening within a fraction of a laser 
period. After the SPW has reached its peak, it is damped 
rapidly. Due to its shorter wavelength, this SPW is even 
more sensitive to thermal damping compared to the SPW 
at the laser fundamental frequency. 

From these results, we deduce that SPW are partic- 
ularly relevant for laser pulses with the duration in a 
certain range. If the pulse is too short, the SPW has 
no time to grow. If the pulse is too long, electron tem- 
perature increases and SPWs are damped rapidly. The 
exact duration depends on the laser and density param- 
eters. In the cases examined here, it was on the range of 
a few up to a few tens laser periods. Note that this is a 
duration range, in which many of today's high intensity 
laser systems work. If SPWs are present, they lead to 
the production of identifiable angular sidebands in the 
harmonics radiation. The angles of the sidebands can 
be estimated by Eq. ([2|. As the /c^-components of the 
sidebands are constant, the angles are greater for lower 
harmonic orders in the optical range and smaller for high 
orders. 
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